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AMPLITUDE DISPERSION AND STABILITY OF 
DISSIPATIVE WEAKLY NON-LINEAR WAVES 

Ferdinand F. Cap 
Theoretical Studies Branch 
Goddard Space Flight Center 
Greenbelt, Maryland 20771 

ABSTRACT 

A perturbation method to solve non-linear dissipative wave equations 
using ideas of Krylov and Bogolyubov is presented. The method is compared 
to Whitham's theory. Dispersion relations for non-linear dissipative 
vaves, including amplitude dispersion, are discussed, Furthermore, 
stability problems of such waves are investigated. 

INTRODUCTION 

Mathematical progress during the last years makes it possible to 
investigate the dispersion of non-linear \<rave equations. Non-linear 
oscillators and waves are important models in hydrodynamics and in 
plasma physics. By the term "non-linear" we want to express that the 
ordinary differential equation describing the propagation of the waves 
are assumed to be non-linear equations. The terms "dispersion" and 
"weakly non-linear" will be discussed later. 

H. Lashingsky^ has discussed in detail the motivation to investi- 
gate mathematical models for non-linear modes in plasmas and he also 

2 

presented such a model starting from a non- linear oscillator equation 

X + (jo 2 x = - eF(x, x) (1) 

To discuss the solutions of this equation, Lashinsky used the method 

■ 

• of averaging in conjunction with the technique of variation of para- 

3 

meters (Krylov-Bogolyubov method ). For e=0. Equation (1) has the 


solution 


x = A cos(a)t + iJj) (2) 

where now A+A(t), for e=j=0. Under the assumptions e<<l, A/A<<u), 

• 

\p/ip«u), which express a weak non- 1 inear i ty , by the so-called method 
of averaging (over one period) the leading terms in a Fourier expansion 
can be found and finally 

• 2tt 

< A > = - 2 ^ / F(A oos0 - a, A sin0)sin0 d0 (3) 

o 

where 0= cut + ip 

The right hand sides are functions of the amplitude A only since 

the time was averaged out over one period. From (3), the amplitude A(t) 

. . 2 

as a function of time t can be computed . Lashinsky also considered 
non-linear waves in bounded plasmas. He investigated the equation 


c 2 


v 2 <p = - CF(<P, A) 


by using expansions of the type 


(4) 


<p = Z a x (t) cj> x (x) (5) 

A 

The orthogonality conditions of the <J>^ and the period averaging 
process result in similar equations like (3) . 

i 

. . *• 4 

The averaging over the period was also used by Tam . He considers 

the propagation of non-linear dispersive waves in a cold non-dissipative 

plasma. He uses a perturbation technique and introduced fast and slow 

variables giving the periodic and the non-linear (averaged) part of the 

5 

solution. Also Luke uses this method in his investigation of the non- 


linear wave equation 


- 3 - 


— <f tt - V 2 4> = F(<(>) 


( 6 ) 


He demonstrates that his perturbation technique is equivalent to 
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Whitham's averaged Lagrangian method. ' ' 

In this paper we extend Whitham's theory of the averaged Lagrangian 
to dissipative waves. Tam 16 points out that the effect of collisional 
dissipation on non-linear dispersion is important. Collisional damping 
may reduce the effect of non-linear instability. 


DISPERSION RELATION AND AVERAGED LAGRANGIAN 

We consider waves <}>(x, t) which satisfy a non-linear partial differ- 
ential equation of such a form that the linear wave ecjuation can be split 
off. Rotating our coordinate system so that the wave vector k points 
into the direction of the x-cixis we write for our wave equation 

— *tt " *xx + ^ + = ' v ' w + N(< M + 

c 2 

* t G(« (7) 

where V' , N and G are non-linear functions, V' = c, b and g are 
constants which may depend on w. We now define a phase surface 

0(x, t) = const (8) 


which has the property that all points (x, t) on it have the same value 
of the wave function <f). From (8) we have 


dO = 0 d + 0 dt = 0 
xx t 


(9) 


so that points moving with the speed 


dx 

it 



( 10 ) 
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see a constant phase 0. Defining wave number and frequency by 

0 X = k, -0 t = a) .. (11) 

we see that (10) is the phase speed. In the three dimensional case 
we have V0 = k, and therefore 

curl k = 0 (12) 


which indicates that wave crests are neither vanishing nor splitting into 

9 

two or more crests. From (11) and (12) the conservation equation of 
wave crests , Equation (13) , follows: 

+ Voj = 0 (13) 

9 10 

In a similar way it can be shown ' that a point moving with the 
group velocity 


<£> 


dco 

dk 


(14) 


sees oo unchanged. After this disgression, we return to Equation (7) . 
Its Lagrangian reads 


L(4>, <P r <P.) = 


x 


1 1 . 2 1 2 b 2 

T 2 ♦t - 2 *x ■ 2 * v(<|,) 


(15) 


Since the systen is dissipative, the Euler Lagrange equations read 


11 , 12 


9 9L , 3 9L 9L _ x , XT/A x , x r \ 

9x 34> 3t 9<f> 9<j) “ g(p t N( V G(cp) (16) 

X L 

where N (<*> t ) may be of the form 

N(* t ) = - e ^ 2 " -1 


( 17 ) 
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where n = 1 1/2, 2, 2 1/2 ... Substituting the Lagrangian (15) into (16) 
we immediately obtain our wave equation (7). We will now make a classifi- 
cation of our wave equation: 

I. N = 0, V' = 0: The equation is linear; k and go are independent of 

x and t. 

1. b = 0, g = 0: No dispersion, no dissipation. 

_ ikx-icot . , /ir> . 

<J) = Ae , go = + ck (18) 


2. b =)= 0, g = 0: Frequency dispersion, no dissipation. 


, _ , ikx-icot , } 

4>-Ae , go - + c / ^ 2^2 


(19) 


Any function co(k) - with exception of the definition (18), is called 
a dispersion relation of a linear wave equation . Such dispersion relations 

may also be found by a Laplace-Four ier transformation of the linear wave 

. . 13 

equation . 

3 . b = 0, g 4 0: Dissipative Case. 


<t> 


Ae 


ikx-icot 

f 



+ c r 


k 2 (t>- g 2 c 2 

4 


( 20 ) 


4. b ^ 0, g 4 0: Dissipative case r 4> like in (20) 


G (go, k) = 


= GO 


- k 2 - b + gico = 0 


( 21 ) 


II. If neither N, G nor V vanish, the wave equation is non-linear and 
solutions cannot be given (except in special cases'^) . If, however, the 


non-linearity is weak , which we define by 

2ir 


2rr 


< G xx > = h ~ f " e xx d0 = 0 ’ <e tt > = h f 0 tt d0 


= 0 


( 22 ) 
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so that 0 = k and 0 = - w become constant "in the average," i.e. when 

X t 2 tt 

appearing under an integral /. . . dO, then a dispersion relation of a 

o 

non-linear wave equation can be derived. It turns out, that is not only 
a function of k (frequency dispersion) but also of the amplitude A (amp- 
litude dispersion) . 

6 7 15 

In the dissipationless case (g = G = N = 0) , Whitham suggested ' ' 
the existence of a variational principle 


3 3X 
3t 30 


+ 


3_ 3*_ 

3x 30 

x 


= 0 


(23) 


for the averaged Lagrangian 

, 2tt 

X- - / L (0, 0 x , © , A)d0 =X(0 x , 0 , A) (24) 

o 

and he assumed that the variation of £ with respect to the amplitude A 
or to the energy E = A 2 /2, i.e. 

X E = 0 or X_ A = 0 (25) 

gives the dispersion relation of the non-linear equation. Equations 

(13), (23) and (25) determine the 3 functions k(x, t) u)(x, t) and A(x, t) 

for exactly linear non-dissipative vibrations and waves the 

15 

condition (25) becomes 


£ = 0 (26) 

This is trivial since L = T-U and the virial theorem states for per- 
iodic motion that <T>=<U>so that the dispersion relation for linear 
motions, i.e. Equation (26) follows from the virial theorem. 


- 7 - 


In the dissipative case a principle analogous to the Whitham principle 
(23) may be found by averaging (16) . A proof of the result will be given 
in the next chapter. We first rewrite (16) using 

<P = <fr(0) (27) 


Substituting into (16) we get 


9 ,d9 9L . 9 ,d0 9L _ 9L d0 

9x l dcf> 9^y 3t ' l d<J> 30^/ 30 d<f> 


d(J) 

*S§ 


+ N - 


dcj> 

d0 


ojG 


(28) 


Averaging now Equation (28) with ~ /. . . d<J> and using (24) gives 


i_ - 2“ *h + 

9x 90 9t 90. 2tt J de v 

X to 


1_ 

2tt 


2tt 2tt , 

/ Nd* - ^ / Gg| d* 


(29) 


0 0 

This is the variational principle which replaces (23) In the dissipative 
case. (13) remains unchanged and an equation replacing (25) will be 
derived in the next chapter in the course of proving (29) . 


PERTURBATION APPROACH 

As mentioned earlier all equations necessary to determine k(x, t) , 

4 5 

w(x, t) and A (x, t) of (18) can be derived ' using a perturbation 
technique and without using a Lagrangian. We now extend this pertur- 
bation technique to the dissipative case. We prefer this technique because 
Whitham' s method gives wrong results in special cases. If one considers 
e.g. (7) , (15) for V=N=G=g= 0, then (25) or (26) do not give the 
right dispersion relation which is given by (19). Whitham thinks ‘ that 


- 8 - 




this restriction of the form of the Lagrangian presumably corresponds 
to the assumption of separability if the Hamilton- Jacobi equation in 
the classical theory of adiabatic invariants. Following Luke we intro- 
duce stretched variables by 

4 

X = ex, T = et (30) 

where e is a small parameter which may give a measure of the wave 

4 

amplitude. In order to include relatively fast local oscillations 
(through the dependence on the variable 0) and to take care of the slow 
variations of A, k and to (through the dependence on the stretched variables 
X a' T) we make the expansions 


<Mx, t) = U(0, X, T) + eUj^G, X, T) 

(31) 

V' (<)>) = V' (U) + eUjV" (U) 

(32) 

Ntt> t ) = N(-U 0 oj) +e(U t -U 10 (o)N' (-U 0 u) 

(33) 

G (<J>) = G(U) + EU G' (U) 

(34) 

Substituting (31) to (34) into (7) we obtain by equating 
powers of £ (and neglecting £ 2 ...) 

the various 

- k 2 )+ bu - go, U g = - V’ (U) + 


N(-UgO)) - G(U)UgU) 

(35) 

which replaces (25) and 


U ]0O & - k2) + “l - 9^10 + V (U) + 


U 10 oN' + UgUt^G' + U lg uG « F 

(36) 


where 
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F = h V + ^GX* + P°e^ + Vx " gU T + 

UrpN' + G (U)U T (37) 

The dispersion relation is contained in (35) . In order to show 
that (36) is equivalent to (29) we first show that = U 0 is a solu- 
tion of the homogeneous equation (37). This is done by derivating 
(35) with respect to 0 and showing that the result is identical to 
(37) for F = 0. In order to solve (36) we then write 


U x = W(0)U 0 (38) 

substituting into (36) one receives 

<|£ - k2 > r w 00 u 0 + “ g^0 u 0 + + 

GW 0 U 0 = F (39) 


Using the identity 


TG ^ U 0 2w 0) “ W 00 U 0 2 + 2VJ Q {J 0e V Q (40) 

and multiplying (39) by U 0 we obtain after elimination of w by (38) and 
integrating 

(£ - k '> < U 10 U 0 - U 1 U 00> = w f < U 10 U 0 - U 1 U 00) (g - N 1 - G)d0 

+ / FU 0 dg (41) 

In order to avoid 17 secular terms proportional to integer powers of 0 
U x and must be bounded. But or are not periodic and are unbounded 
unless the integrals in (41) are bounded for large 0 . Now, if U(and 11 , 
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U T = ~^4Jg) are periodic, then the integrals are bounded. For weak non- 
linearity and weak linear damping we may assume that U is periodic. 
Since any nonperiod’c function multiplied by a periodic function is itself 
a periodic function and since integrals over periodic functions over one 
period are bounded we have from (41) the condition 

/ FU 0 d 0 = 0 (42) 

or using (37) 

If f U 0 f30) + lx (k/ u 0 2 d0) = 

0 0 

2tt 2n 

f (gw - wG)IJ 0 dU + / NdU (43) 

o o 

From (15) and (24) we see that this is equivalent to (29). So now 
(13), (35) - instead of (25) - and (43) determine w(x, t) , k(x, t),A(x, t) . 


APPLICATIONS 

In some special cases the non-linear equations can be solved exactly. In 

5 

these cases 4> (0) is known and (25) can be used. So Luke gives the disper- 
sion relation for the case N=G=g=b=0. This dispersion relation 
is derived from (25) and also frcm (35). Another case which can be integrated 
exactly is b = V' = G = 0. Substituting U 0 = W in (35) and integration 
yields 



k ) 


. 2t r _ 

1 r dW 

2tt 0 g^ 7 + N 


0 

0 


(44) 


For bounciness we have 0-0=1 and (44) gives a dispersion relation. 
If, however, <J>(G) is unknown , ^ cannot be calculated and (29) cannot be used. 
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In these cases either an expansion of can be used, e.g. for stability 
18 

problams or the ordinary differential equations (35) and (36) have 

3 17 

to be solved, e.g. by the Krylov- Bogolyubov method ' or by the Lie 
series method. ^ 

I. Nonlinear conservative equation (N = g = G = 0 

Fran (7) and (27) , or from (35) , we have 

U 00 ( ?F " k2 > + tU + V' (U) = 0 (45) 

An exact solution (by multiplication by Uq ) gives from (24) , (25) 
the dispersion relation (amplitude dispersion) 

- k 2 f dU 2rr = 0 (46) 

° ^2 (E - V) - bU 2 

where E is the energy integration constant. So A = /2E = const. The 
functions u>(x, t) and k(x, t) can be obtained from (13) and (43). The 
result A = const for conservative equations may also be derived by the 
Krylov- Bogol yubov method. For V' = 0, (45) has the generating solution 

A sin (a + ^) , U© = Aaoos(a9 + \ p) (47) 

where A, a, ^ are constants and 



Now, according to Krylov-Bogolyubov we let A-* A(0), (9) and we 

then have from (47) 

Aacos (aGH-q) = As in (aGHip) + A(a+ij)cos(ct0t^) 


(49) 
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• dA 

and therefore (A = ^-) 


Asin (otG+i/j) + A4 >cos (aO+iJ;) = 0 


(50) 


Substituting Uqq = Aotcos(aOti|;) - Aoi(anj>) sin(aOfijj) into (45) gives 


A = - -£ V'oos (aQu/j) 

^ = v,sin ( a0+ ^) 


(51) 

(52) 


Extending V' (U) sin(aGt^) resp V' (U) cos(aO+^) into Fourier series and 
i 211 

integrating J ...dip gives 


2lT 

< A > = ^ — ■=- / V' (U)cos 4d\p 
b 2TT 0 


(53) 


which might be oonpared to (3) and 


2tt 

= /v ' (U)sinW 

Ab 2^r o 


(54) 


Using (47), (53) can be written ' V' (AsiniJ;) dsinij; = C, and we then 
have THEOREM I : For any non-linear conservative wave equation of the 
type (7) the an.olitude A (and therefore the er'gy E) is constant. 
Defining an effective frequency 


ft = a + ip 


(55) 


we may write (47) in the form 


U - A (0) sin (/ft d© +^ q ) 


(56) 


where ib = const. Equation (54) and 

o 


0(x, t) = k(x, t)x - a)(x, t)t 


(57) 
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then give a dispersion relation for see however (46) . 

II. Non-linear Dissipative Case (V' = 0) 

From (35) we now have 

2 

U 00 (^- k 2 ) + bu- g„u = N(-U 0 u) - G(U)U 0 co (58) 

so that V' in (53) and (54) has to be replaced by 

V' ^-gaA(0)acos(a0+^)+ojG(Asin[aGH-i|;] )+N (-uAacos [aGH-^] ) 

Then from (54) we have < \p > =0 and the THEOREM II : For any non-linear 

dissipative wave equation of the of the type (7) the frequency co is 
not modified by the dissipation terns. 

III. Other Problems 

The case b = 0 presents problems here, since in this case the Krylov- 
Bogolyubov method cannot be applied to (35) . Another interesting problem is 

” y t + y xx = by + Xf M (59) 

This equation is of interest in turbulence theory. It can be treated by 
the methods discussed here. 

STABILITY 

In many cases the equation (35) cannot be solved exactly or the wave 

18 

equation is not of the type (7) . In these cases Lighthill proposed for 
conservative systems an expansion of .!£. We would like to follow up a similar 
way for the dissipative case. In general, £ = -i(u), k, A) . We are, however, 
able to eliminate, e.g. , A from one of the three equations determining oj, k, 

A or from the amplitude dispersion relation. With = it (go, k) we rewrite 


(29) by exchanging the order of differentiation 


0 £ - 
tt ox*) 


20 + 
xt oik 


°xx^kk 2lJQ 


( 60 ) 


This is a partial differential equation forO(x, t) . We now expand 
£. Since <)£ . = 0 according to (36) , X. measures exactly the derivation from 

the dispersion relation of the linear equation. In order to make this more 
explicit we expand 


T 


«£(<*>, k) = £(x, k) = ax 2 = a (a)— f [k] f 


(61) 


where (*>= f (k) is the dispersion relation of the linear equation, e.g. (21). 

3 9 3 9 3 

We then have the correspondence - f '-^r, yp where t and k are 

now two independent variables. Since then £ = 0,t = 2ax,/ = 2a, (60) takes 

iC T T 

the form 


®tt + 2f,0 xt + (f,2 ‘ f " T)0 xx 0/3 (62) 

Since f = f (k) , 0 = k(x, t) . This is a quasilinear partial differential 
equation for 0(x, t) . Applying usual character ics method, see e.g. Ref. 20, 
one may write (62) in the form 


AO + B0 . 

XX xt 


+ CO 


tt 


TTQ/a 


(63) 


Then the characteristics are: 

Real - (hyperbolic equation) , if B 2 - 4 AC E 4f "x > 0 
Complex (elliptic equation) , if B 2 - 4 AC = 4f "x < 0 

T^e see that the dissipative te- m Q does not enter into this condition. 
So we have the following THBORJM III: 




■I 







- 15 - 

a. The effect of non-linear terms on the stability behavior is described 
by f " (k) . ( u>- f[k]). The stability behavior of the linear equation, described 
by a) = f (k) is not altered by non-linear terms, if f " (k) . ( u>- f [k] ) > 0. 

If however f" (k) . (u) - f [k] ) < 0, then the non-linear terms may destabilize 
an c herwise stable solution of a linear equation. 

b. The inclusion of dissipative term Q does not by itself modify the 
character of the stability behavior, but the time behavior of unstable and 
stable modes is modified. 

Examples and applications of this theorm will be given in forthcoming 

papers. 
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